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As one step towards a systematic modeling of the electromagnetic (EM) emission from an inspiralling black 
hole binary we consider a simple scenario in which the binary moves in a uniform magnetic field anchored to 
a distant circumbinary disc. We study this system by solving the Einstein-Maxwell equations in which the EM 
fields are chosen with strengths consistent with the values expected astrophysically and treated as test-fields. Our 
initial data consists of a series of binaries with spins aligned or anti-aligned with the orbital angular momentum 
and we study the dependence of gravitational and EM signals with different spin configurations. Overall we find 
that the EM radiation in the lowest £ = 2, m = 2 multipole accurately reflects the gravitational one, with identi- 
cal phase evolutions and amplitudes that differ only by a scaling factor. This is no longer true when considering 
higher I modes, for which the amplitude evolution of the scaled EM emission is slightly larger, while the phase 
evolutions continue to agree. We also compute the efficiency of the energy emission in EM waves and find 
that it scales quadratically with the total spin and is given by E^/M ~ 10~ 15 (M/10 8 M ) 2 (B/10 4 G) 2 , 
hence 13 orders of magnitude smaller than the gravitational energy for realistic magnetic fields. Although 
large in absolute terms, the corresponding luminosity is much smaller than the accretion luminosity if the 
system is accreting at near the Eddington rate. Most importantly, this EM emission is at frequencies of 
~ 10 _4 (10 8 M Q /M) Hz, which are well outside those accessible to astronomical radio observations. As a 
result, it is unlikely that the EM emission discussed here can be detected directly and simultaneously with the 
gravitational-wave one. However, indirect processes, driven by changes in the EM fields behavior could yield 
observable events. In particular we argue that if the accretion rate of the circumbinary disc is small and suffi- 
ciently stable over the timescale of the final inspiral, then the EM emission may be observable indirectly as it 
will alter the accretion rate through the magnetic torques exerted by the distorted magnetic field lines. 



PACS numbers: 

I. INTRODUCTION 

Gravitational-wave (GW) astronomy promises to revolu- 
tionize our understanding of a number of astrophysical sys- 
tems. Several Earth-based detectors (LIGO, Virgo, GEO) are 
already operating at their designed sensitivities and will be 
further upgraded in the coming years. Additionally, space- 
borne detectors are being considered and might become a re- 
ality in the coming decade. The ability to harness the informa- 
tion carried by GWs will soon provide a completely new way 
to observe the universe around us. These detectors, along with 
increasingly sensitive electromagnetic (EM) telescopes, will 
provide insights likely to affect profoundly our understanding 
of fundamental physics and the cosmos (see e.g. fl] EJ for a 
recent and detailed discussions of the astrophysics and cos- 
mology that will be possible with the detection of GWs) 

Among the most promising sources of detectable GWs are 
systems composed of binary black holes which, as they come 
together and merge, radiate copious amounts of energy in the 
form of gravitational radiation. When these black holes are 
supermassive, i.e., with masses M > 10 6 Mq, the cosmolog- 
ical and astrophysical conditions leading to their formation 
will be such that prior to the merger they will be surrounded 
by a gas or plasma that could also radiate electromagnetic ally. 
Indeed, within the context of galaxy mergers, such a scenario 
will typically arise as the central black hole in each of the col- 



liding galaxies sinks towards the gravitational center, eventu- 
ally forming a binary. In such a process, the binary will gener- 
ically find itself inside a circumbinary disc which might be a 
catalyst of observable emissions as it interacts with the black 
holes [3]|U 0]. Within this context, several possibilities are ac- 
tively being investigated. Among these, several studies have 
(with varying degrees of approximations) concentrated on un- 
derstanding emissions by the disc due to the interaction with a 
recoiling, or a mass-reduced final black hole [6 7 8 9, iOlfTTIl 
and remnant gas around the black holes lfl2l [T3l [141 [T5l . 

A further intriguing possibility for such interaction is 
through EM fields and in particular through the gravitomag- 
netic deformation of magnetic fields, anchored to the disc, 
around the central region where the black holes inspiral and 
merge. As the black holes proceed in an ever shrinking orbit 
towards their ultimate coalescence, they will twist and stir the 
EM fields and thus affect their topology. Moreover, the space- 
time dynamics might impact the fields in such a strong way 
as to generate EM energy fluxes which may reach and impact 
the disc and/or affect possible gas in the black holes' vicin- 
ity. Last, but certainly not least, as the merger takes place the 
system could acquire (at least at much later times) a configu- 
ration where emissions through the Blandford-Znajek mecha- 
nism may take place [ 16 1 

The initial step towards understanding this system was 
taken in ifTTl [18 1 and highlighted the possible phenomenol- 
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ogy in the system. Although this first investigation was re- 
stricted to the simplest case of an equal-mass, nonspinning 
binary, it provided a proof of principle examination of possi- 
ble EM counterparts resulting from binary mergers. The work 
presented here extends the study carried out in ifTTl [181 on 
several fronts. First, it considers an broader class of binaries, 
with and without net spins which are aligned/counter-aligned 
with the orbital angular momenta (in the case of supermassive 
black holes, these configurations are indeed the most likely to 
be produced 1 19 20]). Second, it studies binaries from larger 
separations, thus allowing for a clearer modeling of the inspi- 
ral. Third, it examines more closely the energy emissions in 
both channels and compares the total output's dependence on 
final black hole spin. Finally, it assesses the direct detectabil- 
ity of the EM radiation produced, and considers the astrophys- 
ical impact it may have on the surrounding accretion disc. In 
addition, while sharing the same astrophysical scenario, the 
methodological approach adopted here also differs from that 
of |[T7l[T8l in at least three important ways. First, we here ex- 
plicitly impose the test-field limit by setting to zero the stress- 
energy tensor on the right-hand-side of the Einstein equations. 
Second, we employ a different formulation of Einstein equa- 
tions, i.e., a conformal transverse-traceless one in place of a 
generalized-harmonic one. Third, we use a distinct compu- 
tational infrastructure based on Cactus and Carpet rather 
than the HAD code which was instead used in IfTTl [TEl ). The 
excellent agreement in the phenomenology observed, in cases 
considered in both works, gives a strong further indication of 
the correctness of the results obtained. 

Our analysis shows that the EM-wave emission in the low- 
est £ = 2, to = 2 multipole reflects in a spectacular way the 
gravitational one and that the phase and amplitude evolution 
differ only by a scaling factor. This is no longer true when 
considering higher I modes, for which the amplitude evolu- 
tion of the EM emission is slightly larger. We also find that the 
efficiency of the energy emission in EM waves scales quadrat- 
ically with the total spin and is given, for realistic magnetic 
fields, by E™£/M ~ 10" 15 (M/10 8 M q ) 2 (B/10 4 G) 2 . 
For expected fields of 10 4 G this energy loss is 13 orders 
of magnitude smaller than the gravitational one. In addition, 
the corresponding EM luminosity is much smaller than that 
of accretion, if the system is accreting at near the Edding- 
ton rate. Most importantly, this emission is at frequencies 
of - 1CT 4 (1O 8 M /M) Hz, which lies outside the range of 
astronomical radio observations. As a result, it is highly un- 
likely that the EM emission discussed here can be detected di- 
rectly and simultaneously with the GW one. Other processes 
however could be affected by this flux of EM energy and pro- 
duce detectable effects. For instance, if the accretion rate of 
the circumbinary disc is small and sufficiently stable over the 
timescale of the final inspiral, then the EM emission may be, 
in particular, observable indirectly as it will alter the accre- 
tion rate through the magnetic torques exerted by the distorted 
magnetic field lines. 

This work is organized as follows: Sect. [II] presents an 
overview of both Maxwell and Einstein equations as imple- 
mented within the BSSNOK formulation of the Einstein equa- 
tions. The physical and astrophysical setup adopted in our 



simulations is discussed in Sect. [Ill] while Sect. [W] is dedi- 
cated to the analysis of the EM fields for isolated black holes, 
either nonspinning or with spin aligned with the orbital angu- 
lar momentum. Sect. [V]collects instead our results for the dif- 
ferent binaries considered, while the assessment of the astro- 
physical impact of the EM emission is presented in Sect. [VT] 
We conclude in Sect. IVIII with final comments and discus- 
sions. 



II. THE EVOLUTION EQUATIONS 

We solve the Einstein-Maxwell system to model the inter- 
action of an inspiralling black-hole binary with an externally 
sourced magnetic field in an electro-vacuum spacetime. More 
specifically, we solve the Einstein equations 



(1) 



where g^ v , R^ v and are the metric, the Ricci and the 
stress-energy tensor, respectively, together with an extension 
of Maxwell equations in absence of currents ET1I221 . 



(2) 
(3) 



when written as conservation laws for the Faraday tensor F^ v 
and of its dual *F )JjV . Note that we have here introduced two 
extra scalar fields W and </>, which are initially zero and whose 
evolution drives the system to a satisfaction of the EM con- 
straints (see discussion below). The two systems are coupled 
through the stress-energy tensor, 
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Note also that since our stress-energy tensor is not identi- 
cally zero, the binary is not in vacuum, at least within a strict 
general-relativistic sense. However, because the sources of the 
EM fields are not part of our description and we simply con- 
sider EM fields, our setup will be that of an "electromagnetic - 
vacuum", which will simply refer to as "vacuum". 



A. The Einstein Equations 

The numerical solution of the Einstein equations has been 
performed using a three-dimensional finite-differencing code 
solving a conformal-traceless "3+1" BSSNOK formulation of 
the Einstein equations (see [23] for the full expressions in vac- 
uum and ll24l for spacetimes with matter) using the Cactus 
Computational Toolkit 11251 and the Carpet ll26l adaptive 
mesh-refinement driver. Recent developments, such as the use 
of 8th-order finite-difference operators or the adoption of a 
multiblock structure to extend the size of the wave zone have 
been recently presented in l23l l27l . Here, however, to limit 
the computational costs and because a very high accuracy in 
the waveforms is not needed, the multiblock structure was not 
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used. For compactness we will not report here the details of 
the formulation of the Einstein equations solved or the form 
of the gauge conditions adopted. All of these aspects are dis- 
cussed in detail in ll23l . to which we refer the interested reader. 

B. The Maxwell Equations 

Maxwell equations |2])-([3]) take a more familiar form when 
represented in terms of the standard electric and magnetic 
fields. These are defined by the following decomposition of 
the Faraday tensor 

F ?u = t^E" -t v E» + e^ afi B a t p , (5) 
*F» V = t»B v -t v B* -e^"! 3 Eatp ; (6) 

where is the unit time vector associated with a generic nor- 
mal observer to the hypersurfaces. The vectors E^ and B^ L are 
the (purely spatial, E^t^ = B^t^ = 0) electric and magnetic 
fields measured by such observer. 

As mentioned above, we adopt an extended version of 
Maxwell's equations which introduces two extra scalar fields 

and <1>. This extension induces evolution equations for the 



EM constraints (ViE % = = V,-B l ) described by damped 
wave equations and so control dynamically these constraints. 
In terms of E 11 and B^ L the 3+1 version of |2]|3]l results, 

V t E l - e ijk V.j (aB c ) + aj ij Vj^ = aKE l , (7) 

V t B l + e ijk Vj (aE c ) + aY 1 V J ^ = aK B l , (8) 

D(f + aV i £ , = -aK$, (9) 

V t ^ + aV i B i = -an®. (10) 

where V t = (dt — Cp) and Cp is the Lie derivative along the 
shift vector /3. Exploiting that the covariant derivative in the 
second term of (|8]-|9|l reduces to a partial one 

e ijk VjB k = e ijk (dj + T c k0 B k ) = e ijk djB k , (11) 

and using a standard conformal decomposition of the spatial 
3 -metric 

7*3 = e 40 7ij ; = (12) 
one obtains the final expressions 



V t E l - e l > k [ (dj a)% k B c +a(A % k 0, 4> + dj % k ) B c + a % k d 3 B c }+a e~ A * f j V^^aKE 1 , (13) 

V t B l + e^ k [ (dj a ) % k E c + a ( 4 7cfc dj <j> + dj % k ) E c + a % k dj E c ]+a f •> V j $ = aKB i , (14) 

Di* + aV,£ l = -«K*, (15) 

V t $ + a ViB % = -an <$> . (16) 



Notice that the standard Maxwell equations in a curved back- 
ground are recovered for W = <f> = 0. The W and <f> scalars 
can then be considered as the normal-time integrals of the 
standard divergence constraints 



ViE 1 = , V,B l = 



(17) 



These constraints propagate with light speed and are damped 
during the evolution. 

As mentioned above, the coupling between the Einstein and 
Maxwell equations takes place via the inclusion of a nonzero 
stress-energy tensor (cf. the set of equations presented in [24 |) 
for the EM fields and which is built in terms of the Faraday 
tensor as described in Q. More specifically, the relevant com- 
ponents of the stress-energy tensor can be obtained in terms of 
the electric and magnetic fields, that is 
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-EiEj — BiBj + — 7y (E 2 + B 2 ) 



(18) 
(19) 



where E 2 = E k E k and B 2 = B k B k . The scalar compo- 
nent t can be identified with the energy density of the EM 



field (i.e., p ADM in l24l ) and the energy flux Si is the Poynt- 
ing vector. As already discussed in the Introduction, we stress 
again that the EM energies considered here are so small when 
compared with the gravitational ones that the contributions of 
the stress-energy tensor to the right-hand-side of the Einstein 
equations (fT} are comparatively negligible and thus effectively 
set to zero 



C. Analysis of Radiated Quantities 

The calculation of the EM and gravitational radiation gen- 
erated during the inspiral, merger and ringdown is arguably 
the most important aspect of this work as it allows us to es- 
tablish on a firm basis whether a clear correlation should be 
expected between the two forms of radiation. We compute 



1 The fully coupled set of the Einstein-Maxwell equations was considered 
in II 711 1 81 and the comparison with the results obtained here suggests that 
for the fields below 10 8 G, the use of the test-field approximation is justi- 
fied. 
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the gravitational radiation via the Newman-Penrose curvature 
scalar ^4 defined as 

* 4 = -C^s^m^n^fh 5 , (20) 

where C a ^s is the Weyl curvature tensor, which is projected 
onto a null frame, {l,n,m,fh}. In practice, we define an 
orthonormal basis in the three space (r, 0,(p), centered on the 
Cartesian grid center and oriented with poles along z. The 
normal to the slice defines a time-like vector t, from which 
we construct the null frame 



1 



V2 



(t-r), 



•=^(* + r), m = -J=(0-i<£). 

(21) 

We then calculate ^4 via a reformulation of ( |20| in terms of 
ADM variables on the slice I 



where 



Ri 



*4 = a 



KKa + Ki k K k 



kl 



We note that the we have also implemented an indepen- 
dent method to compute gravitational radiation via the mea- 
surements of the nonspherical gauge-invariant metric pertur- 
bations of a Schwarzschild black hole J29] |30l [31] (A re- 
view on the basic formalism can be found in 11321 and 
refs. Il29l l30l I3TI provide examples of the applications of 
this method to Cartesian-coordinate grids). For compactness, 
hereafter we will limit our discussion to the gravitational ra- 
diation measured in terms of the curvature scalar 

In a similar manner, the EM radiation can be measured us- 
ing the Newman-Penrose scalar $2 defined as 

$ 2 = F^fh^n, , (24) 

with the same same tetrad used for ^4, allowing to measure 
outgoing EM radiation. (Possible gauge effects, as those dis- 
cussed in 11331 . have been seen to be negligible in 1181 . and the 
results here agree with those). 

Using the curvature scalars ^4 and $2 it is also possible to 
compute the energy carried off by outgoing waves at infinity. 
More specifically, the total energy flux per unit solid angle can 
be computed directly as 



(22) 



(23) 
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(25) 



(26) 



dtdE 

Representing now ^4 and $2 via a standard decomposition 
into spherical harmonics using the spin-weights —2 for ^4 
and —1 for $2, 
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, we obtain the final expressions 
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where A lrn and B lm are the coefficients of the spherical har- 
monic decomposition of ^4 and $2, respectively. 



III. PHYSICAL AND ASTROPHYSICAL SETUP 

As mentioned in the introduction, the astrophysical sce- 
nario we have in mind is motivated by the merger of super- 
massive black holes binaries resulting from galaxy mergers. 
More specifically, we consider the astrophysical conditions 
that would follow the merger of two supermassive black holes, 
each of which is surrounded by an accretion disc. As the 
merger between the two galaxies takes place and the black 
holes become close, a "circumbinary" accretion forms and 
reaches a stationary accretion phase. During this phase, the bi- 
nary evolves on the timescale of the emission of gravitational 
radiation and its separation progressively decreases as gravita- 
tional waves carry away energy and angular momentum from 
the system. This radiation-reaction timescale is much longer 
than the (disc) accretion timescale, which is regulated by the 
ability of the disc to transport outwards its angular momen- 
tum (either via viscous shear or magnetically-mediated insta- 
bilities). As a consequence, for most of the evolution the disc 
slowly follows the binary as its orbit shrinks. However, as the 
binary separation becomes of the order of ~ 10 5 — 10 6 M , the 
radiation-reaction timescale reduces considerably and can be- 
come smaller than the disc accretion one. When this happens, 
the disc becomes disconnected from the binary, the mass ac- 
cretion rate reduces substantially and the binary performs its 
final orbits in an "interior" region which is essentially devoid 
of gas |3]|4]|5]- This represents the astrophysical setup of our 
simple model. 

We introduce a coupling between the binary and the disc via 
a large-scale magnetic field which we assume to be anchored 
to the disc, whose inner edge is at a distance of ~ 10 3 M and 
is effectively outside of our computational domain, while the 
binary separation is only of ~ 10 M, where M is the total 
gravitational mass of the binary. We note that although the 
large-scale magnetic field is poloidal, it will appear as essen- 
tially uniform within the "interior region" where the binary 
evolves and which we model here. As a result, the initially 
magnetic field adopted has Cartesian components given sim- 
ply by B i = (0, 0, B ) with B Q M = 10~ 4 in geometric units 
or B a - 10 s G for a binary with total mass M = 10 8 M Q . 
Furthermore, because we consider an electromagnetic vac- 
uum, the charges, electric currents and the initial electric field 
are all assumed to be zero, i.e., E % = 0. 

We note that although astrophysically large, the initial mag- 
netic field considered here has an associated EM energy which 
is several orders of magnitude smaller than the gravitational- 
field energy. As a result, any effect from the EM field dynam- 
ics on the spacetime itself will be negligible and so the EM 
fields are treated here as test-fields. The case of stronger mag- 
netic fields and their consequent impact on the spacetime will 
be presented in a forthcoming work. 



IV. ISOLATED BLACK HOLES 

We first study isolated black holes, both as a check of our 
implementation and to analyze the interaction of the chosen 
external initial magnetic field with the spacetime curvature 
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FIG. 1 : Recovery of Wald's solution for an isolated Kerr black hole 
with dimensionless spin a = 0.7. Shown in the top panel are the 
values of the electric field as measured at different distances from 
the origin; since Er 2 ~ Bo J, the different lines should overlap at 
late times if the magnetic field is uniform which is evident in the 
figure. Shown in the bottom panel is the ratio of the electric and 
magnetic fields which is proportional to the black hole spin only. 
Note the transitory state until t ~ 70M, when the solution reaches a 
stationary state. 



generated by the black holes. The initial magnetic field in all 
simulations is uniform with strength Bq and aligned with the 
z-axis, while the initial electric field is zero everywhere. Al- 
though this solution satisfies the Maxwell equations trivially, 
it is not a stationary solution of the coupled Einstein-Maxwell 
system for the chosen black hole initial data. The solution 
thus exhibits a transient behavior and evolves towards a time- 
independent state given by a solution first found by Wald (34). 
One important feature of Wald's solution is that in the case 
of spinning black holes, a net charge (and hence a net elec- 
tric field) will develop as a result of "selective accretion" and 
whose asymptotic value is simply given by Q — 2B J. Al- 
though this charge is astrophysically uninteresting, being lim- 
ited to be Q/M < 2B M ~ 1.7 x 1CT 20 B (M/M Q ) G El 
for a Kerr black hole with J/Al 2 < 1, it represents an excel- 
lent testbed for our numerical setup. 

To validate the ability of the code to recover this analytic 
solution we have performed several tests involving either a 
Schwarzschild black hole or Kerr black holes with dimension- 
less spin parameters a = J/M 2 = 0.7 (this value chosen as it 
is close to the final spin values resulting from the merger sim- 
ulations covered in Sect.|IV[). In this latter case, the spin vector 
was chosen to be either parallel to the background magnetic 
field, i.e., with J 1 = (0, 0, J) or orthogonal to it, i.e., with 
J 1 = (J, 0, 0). As expected, the early stages in the evolution 
reveal a transient behavior as the EM fields rearrange them- 
selves and adapt to the curved spacetime reaching a stationary 



configuration after about ~ 70 M, The electric field, in par- 
ticular, goes from being initially zero to being nonzero and 
decaying radially from the black hole. 

Although the original solution found by Wald was ex- 
pressed in Boyer-Lindquist coordinates, there is a simple way 
to validate that our gauge is sufficiently similar (at least at far 
distances) and that the numerical solution approaches Wald's 
one for an isolated black hole in a uniform magnetic field. 
This is shown in Fig. [T] which reports the time evolution of 
the EM fields E and B for a simulation of the Kerr black hole 
with spin a = 0.7 M aligned with the magnetic field. In partic- 
ular, the top panel shows the time evolution of the electric field 
when the latter is rescaled by the radial positions where it is 
measured, i.e., Er 2 with r = 4 M, 8 M , 16 M and 24 M. Be- 
cause of E tx B a J/r 2 , one expects the different lines to be on 
top of each other. This is clearly the case for the data extracted 
at r = 16 M and 24 M, but it ceases to be true for the data at 
r = 4 M, 8 M, for which the magnetic field and gauge struc- 
ture are strongly influenced by the black-hole geometry. In- 
terestingly, however, in this strong-field region near the black 
hole another scaling can be found and which is closely related 
to one expressed by Wald's solution. In particular, the radial 
dependence of the magnetic field can be factored out by con- 
sidering the ratio of the electric and magnetic field which, in 
Wald's solution, should be proportional to the black-hole spin 
only. The bottom panel of Fig.[T]shows therefore the evolution 
of (E/B)r 2 ~ J which is indeed a constant at all the radial 
positions as shown by the good overlap among the different 
curves. We find that this scaling can be used as an effective 
test which is valid at all radial positions. These observations, 
together with the clear approach to a stationary configuration 
indicate the asymptotic (in time) behavior is indeed described 
by Wald's solution. 

In order to obtain a more intuitive picture of the differ- 
ent solutions for isolated black holes, we now turn our atten- 
tion to the structure of the electric and magnetic fields them- 
selves. While those field lines are gauge-dependent, they can 
be used to determine the effect of the spin orientation of the 
black holes on the solution. Fig. [2] shows therefore the three- 
dimensional (3D) EM field configurations at late simulation 
times when the solution has settled to a stationary state for ei- 
ther a Schwarzschild black hole (left panel), or for Kerr black 
holes with spin aligned (central panel) or orthogonal to the 
magnetic field (right panel). Note that in all of the panels, 
the magnetic field lines are bent by the black hole geome- 
try. The appearance of toroidal electric field in the case of a 
nonspinning black hole does not contradict Wald's solution, 
for which it should be identically zero. It is due to the non- 
vanishing radial shift vector which, when coupled with the 
vertical magnetic field, leads to a toroidal magnetic field ifTHl . 
Finally, note that whenever the black hole is rotating, together 
with the gauge-induced toroidal electric field, there appears 
also a poloidal component which is induced by the gravito- 
magnetism (or frame-dragging) of the rotating black hole and 
whose detailed geometry depends on the relative orientation 
of the spin with respect to the background magnetic field. For 
compactness we do not report here the EM field configuration 
for a rotating black hole with spin anti-aligned with respect to 



FIG. 2: Left panel: Magnetic (blue) and electric field (red/magenta) field lines at t = 200 M for a Schwarzschild black hole. Central panel: 
the same as in the left panel but for a Kerr black hole with spin a = J/M 2 = 0.7 aligned with the magnetic field, i.e., J 1 = {0, 0, J}. 
Right panel: the same as in the center panel but for a Kerr BH with spin a = J/M 2 = 0.7 which is orthogonal to the magnetic field, 
i.e., ,T = {J, 0, 0}. Indicated with black surfaces are the apparent horizons. 



the magnetic field. It is sufficient to remark that the solution 
shows the same behavior as the aligned case, with a simple 
reversal in the direction of the spin-induced effects. 

To gain some insight on the influence of the black hole spin 
and orientation on the EM field lines, it is useful to exploit 
the phenomenological description offered by the "membrane 
paradigm" ||35l . In such approach, the horizon of a rotating 
black hole is seen as a one-way membrane with a net a sur- 
face charge distribution which, for the case of aligned spin 
and magnetic field, has negative values around the poles while 
positive ones around the equator. The resulting behavior is 
therefore the one shown in Fig. [3] where the magnetic and 
electric field lines for the Kerr black hole with spin aligned 
with the magnetic field are presented on the y = 0-plane. 
The left panel, in particular, offers a large-scale view of the 
EM fields, which is however magnified on the right panels to 
highlight the behavior of the fields near the apparent horizons. 
Finally, shown in Fig. |4] are the magnetic and electric-field 
lines on the y = 0-plane for the Kerr black hole with spin or- 
thogonal to the the magnetic field, i.e., along the a>axis. Note 
that while the differences in the magnetic field configurations 
in Fig. [5] and [4] are small and difficult to observe even in the 
zoomed-in version of the figures, the differences in the elec- 
tric fields are instead significant and related to the different 
spin orientations. 



V. BINARY BLACK HOLES 

We next extend the considerations made in the previous sec- 
tion to a series of black hole binaries having equal masses and 
spins that are either aligned or anti-aligned with the orbital 
angular momentum. 



A. Initial Data and Grid Setup 

We construct consistent black-hole initial data via the 
"puncture" method as described in ref. (36). We consider 
equal mass binaries with four different spin configurations be- 
longing to the sequences labeled as "r" and "s" along straight 
lines in the (01,02) parameter space, also referred to as the 
"spin diagram" |[37l [38]. These configurations allow us to 
cover the basic combinations for the alignment of the spin of 
the individuals black holes with respect to the magnetic field, 
while keeping the dimensionless spin parameter of the single 
black holes constant among the different binaries considered. 
Furthermore, it allows us to study the impact that the final 
black hole spin has on the late stages of the merger. 

We note that similar sequences have also been considered 
in (39] 00] |4T] g2l mi but have here been recalculated both 
using a higher resolution and with improved initial orbital pa- 
rameters. More specifically, we use post-Newtonian (PN) evo- 
lutions following the scheme outlined in [44|, which provides 
a straightforward prescription for initial-data parameters with 
small initial eccentricity, and which can be interpreted as part 
of the process of matching our numerical calculations to the 
inspiral described by the PN approximations. The free pa- 
rameters of the puncture initial data we fix are: the puncture 
coordinate locations Ci, the puncture bare mass parameters 
m„ the linear momenta p,, and the individual spins Si. The 
initial parameters for all of the binaries considered are col- 
lected in Table[I] The initial separations are fixed at D = 8 M 
with the exception of the s_6 binary having an initial sepa- 
ration of D = 10 M. Here M is the total initial black hole 
mass, chosen as M = 1 (note that the initial ADM mass of 
the spacetime is not exactly 1 due to the binding energy of 
the black holes), while the individual asymptotic initial black 
hole masses are therefore M$ = 1/2. In addition, we choose 
the initial parameters for the EM fields to be B % = (0, 0, Bq) 
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FIG. 3: Left panel: Large-scale magnetic and electric field lines on the plane y = and at t — 200 M for a Kerr black hole with spin 
J/M 2 — 0.7 aligned with the magnetic field, i.e., along the z-axis. Indicated with blue circles are the apparent horizons. Right panel: The 
same as on the left panel but on a smaller scale to highlight the fields structure in the vicinity of the black hole. 
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FIG. 4: Left panel: Large-scale magnetic and electric field lines on the plane y = and at t — 200 M for a Kerr black hole with spin 
J/M 2 = 0.7 orthogonal to the magnetic field, i.e., along the a;-axis. Indicated with blue circles are the apparent horizons. Right panel: The 
same as on the left panel but on a smaller scale to highlight the fields structure in the vicinity of the black hole. 



with B ~ 10~ 4 /A/ - 1O 8 (1O 8 M /M) G and E i = 0. The 
setup for the numerical grids used in the simulations consists 
of 9 levels of mesh refinement with a fine-grid resolution of 
Ax/M — 0.02 together with fourth-order finite differencing. 
The wave-zone grid has a resolution of Ax/M — 0.128 and 



extends from r = 24 M to r = 180 M, in which our wave 
extraction is carried out. The outer (coarsest) grid extends to 
a spatial position which is 819.2 M in each coordinate direc- 
tion. 
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FIG. 5: Electric field lines on the plane y = for the ro, so and S6 configurations at t — 123, 155 and 246 M, respectively. Left panel: 
Large-scale structure of the EM fields around the apparent horizons (blue circles). Right panel: The same as on the left but on a smaller scale 
to highlight the field structure in the vicinity of the black holes. Note that an additional magnification is applied to the black hole "on the right" 
so as to highlight the change of sign in the electric field near the horizon, i.e., at x ~ 3 M. 




FIG. 6: Electric (red/magenta) and magnetic field lines (gray) in 3D for the s& binary during inspiral when both black holes are still far 
separated at time t = 328 M (left panel), and after the merger at t = 690 M (right panel). 



B. Binary Evolution and Spin Dependence 



As mentioned above, we consider configurations where 
both black holes have equal mass and the individual black 
hole spins are either aligned or anti-aligned with the magnetic 
field (and orbital angular momentum). We thus consider a set 
of three different spinning binaries, as well as a nonspinning 
binary, which we take as a reference (cf. Table 



One feature of our simulations, that was already analysed 



for single black holes in Sect. IV and is of even greater inter- 



est for binaries, is the structure of the EM field lines induced 
by the spacetime dynamics around the black holes. The field 
line configurations, in fact, change considerably throughout 
the course of our simulations. When there is a large separation 
between the orbiting black holes, the electric field structure in 
both nonspinning and spinning binary systems is dominated 
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TABLE I: Binary sequences for which numerical simulations have been carried out, with various columns referring to the puncture initial 
location ±x/M, the mass parameters rrii/M, the dimensionless spins a;, the initial momenta and the normalized ADM mass A/ ADM = 
M ADM /M measured at infinity. (See refs. [37, 38 1 for a discussion of the naming convention). 





±x/M 


mi/M 


m 2 /M 


ai 


a 2 


(Px, Py)l = ~(Px, Py)l 


^ADM 


S-6 


5.0000 


0.4000 


0.4000 


-0.600 


-0.600 


(0.001191,-0.100205) 


0.9873 


ro 


4.0000 


0.4000 


0.4000 


-0.600 


0.600 


(0.001860, -0.107537) 


0.9865 


so 


4.0000 


0.4824 


0.4824 


0.000 


0.000 


(0.002088,-0.112349) 


0.9877 


S6 


4.0000 


0.4000 


0.4000 


0.600 


0.600 


(0.001860, -0.107537) 


0.9876 




FIG. 7: Left panel: GWs as computed from the (2, 2) mode of ^4 for the different binaries reported in Table]!] Right panel: The same as in 
the left panel but for the EM waveform as computed from $2. 



by the orbital motion of the individual black holes. In particu- 
lar, an inspection of the electric field vector along a line join- 
ing their centers indicates an outward radial dependence. This 
can be understood from the phenomenological interpretation 
suggested by the membrane paradigm and has been observed 
already in ff8l . Namely, as the black holes move in a direc- 
tion which is essentially orthogonal to the magnetic field, an 
effective quadrupolar charge separation develops on the hori- 
zons with effective positive charges at the poles and negative 
ones on the equator, thus inducing an electric field emanating 
from each black hole. This induced quadrupolar electric field 
is therefore reminiscent of the one produced by a conductor 
moving through a uniform magnetic field as the result of the 
Hall effect. 

It is interesting to note that while the differences in the 
magnetic field lines among the various binaries considered are 
rather small, the differences in the electric fields show signif- 
icant variation across the spin configurations. This is illus- 
trated in Fig. [5] which shows the electric field lines at differ- 
ent scales of interest with respect to the black holes for two 



spinning binary black hole systems and the nonspinning bi- 
nary on the y = 0-plane. Here we choose to concentrate on 
the configurations with the spins up/up (i.e., Sg) and with the 
spin up/down (i.e., ro) since the configuration with the spins 
down/down (i.e., s-e) shows the same field-line structure as 
the up/up case. In particular, the left panel of Fig. |5]reports the 
field-line structure on a scale which is much larger than that 
of the horizons and that clearly shows the quadrupolar nature 
of the field. At the same time, the right panel offers a magni- 
fied view of the same binaries on scales which is comparable 
with those of the horizons. In this way it is possible to find the 



properties of the electric field already discussed in Sect IV for 



isolated black holes also in case of binary black holes. Addi- 
tionally the various spin configurations lead to different small- 
scale properties of the field. More specifically, while the field 
lines of the ro and sq configurations have a similar structure 
even in the magnified plot, the binary with the aligned spins 
s 6 shows a more complex structure in which the electric field 
changes sign near but outside the horizon, namely at x ~ 3 M 
and which corresponds approximately to a distance d ~ 2r AH , 
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t - r (M) t - r (M) 

FIG. 8: Amplitude and phase evolution of the main I = m = 2 modes for the Weyl scalar $2 and the first time integral of ^4 (i.e., $4), 
relative to the ss configuration. The plots show the data in retarded time t — r for a detector located at r = 100M. While the I = m = 2 
modes show the same amplitude (up to a scale factor) and phase evolution, this does not apply to modes with higher I. For the t — 3, m = 2 
and £ = 4, m = 2 modes, the phase evolution is still identical but the amplitude no longer does not differ only by a constant scale factor. 



with r AH the mean radius of the apparent horizon. This addi- 
tional property of the electric field could be related to the loca- 
tion of the ergosphere (which has not been computed in these 
simulations) and may be seen as a response of the electric field 
to the additional charge separation induced on the black hole 
horizon and which leads to a greater distortion and twisting of 
the field lines in this region. 

Although it is not trivial to disentangle how much of this 
behaviour of the electric field depends on the gauges used, 
the complex structure of the electric fields, and which varies 
considerably through the late inspiral and the merger of the 
binary, may lead to interesting dynamics and to the extraction 
of energy via acceleration of particles along open magnetic 
field lines or via magnetic reconnection. To better illustrate 
the complex field structure, Fig.|6]offers 3D views of both the 
electric and magnetic fields for the sq binary. In particular, 
in the left panel of Fig. [6] we show the electric and magnetic 
field lines as well as the apparent horizons when the binary 
is inspiralling (i.e., at t — 328 M) and again observe the su- 
perposition of two effects: the overall orbital motion of the 



black holes causing the large scale structure of the electric 
field lines (highlighted in a magenta color); and the effect of 
the black hole spin (in red), which causes additional dragging 
in the electric field lines close to the apparent horizons. In 
the right panel, on the other hand, we present the late-time 
(i.e., att — 690 M) state of the solution which, as expected, 
agrees well with the field line configurations presented for the 
Kerr black hole with spin aligned with the magnetic field in 
Sect.HV] 

We next switch our attention to discussing how the different 
BH spin configurations affect the emission of EM radiation. 
This requires a careful analysis of the radiative properties of 
the solution in both the EM and gravitational channels. We 
first focus on the two types of waveforms and Fig. [7] illus- 
trates the correlation between the two emissions by showing 
the time-retarded waveform of the principal mode, i.e., the 
I = m = 2 of the spin-weighted spherical harmonic ba- 
sis (note that ^4 and $2 have spin weight —2 and — 1, re- 
spectively), for all different spin configurations to compare 
gravitational and EM waveforms directly. While both EM 
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FIG. 9: The total energy flux per unit solid angle in terms of GW waves (left panel) and of EM waves (right panel); clearly they differ only up 
to a scaling factor. The different lines refer to the different binaries reported in Table [I] 



and GW radiation show the same characteristics in the domi- 
nant mode, we note that small differences arise when compar- 
ing the waveforms of the individual spin configurations more 
carefully with each other in the two channels. Note that the 
waveform for the binary s_6 has a larger number of cycles 
only because it merges very rapidly (the total angular momen- 
tum is smaller because the total spin is anti-parallel to the or- 
bital angular momentum) and thus it has been evolved from a 
larger initial separation D = 10 M; all the other binaries have 
the same initial separation D = 8 M. A closer inspection of 
Fig.|7]reveals that the amplitude evolution of the £ = to = 2- 
mode for the different spin configurations differs when com- 
pared in the two channels. As an example, while in the GW 
channel the amplitude in the £ = 2, to = 2-mode decreases 
when going from the ro -configuration over to the sq and s_6 
configurations, the amplitude remains nearly constant in the 
EM channel. This reveals that there are additional contribu- 
tions in the EM emission coming from the higher-order modes 
(see Fig.|9]and the discussion below) 

To further evaluate the correlation between the EM and the 
gravitational radiation, we now turn our attention to the am- 
plitude and phase evolution of the main contributing spheri- 
cal harmonic modes. Since radiated energy fluxes are given 
by $2 and the time integral of ^ 4 we here compare <& 2 with 
\&4 = f ty^dt'. For briefness we only highlight the results 
obtained for the s 6 configuration, since this shows the high- 
est amount of energy being radiated in both EM and GWs, 
and because our remarks apply also to the other configura- 
tions. Since the main contributions to the radiated energy in 
the EM channel arise from the £ — 2,3,4. m = 2-modes, 
we limit our analysis in this section to those modes only. In 



order to obtain a better understanding of the correlation in the 
radiation coming from the two channels, we analyse the am- 
plitude and phase of the main contributing modes individu- 
ally. Fig. [8] shows the amplitude and phase evolution of the 
£ = 2, 3, to = 2-modes in both channels. Clearly, the 
£ = to = 2 modes show the same phase evolution (cf. the 
left panels) in the two forms of radiation, as expected given 
that the EM emission is essentially driven by the orbital mo- 
tion of the binary. Furthermore, the amplitude evolution in the 
I = 2, to = 2 modes of both emissions are also simply related 
by a constant, time-independent, factor. 

Although a simple scaling factor in the evolution of "J 4 and 
$2 appears for all of the different binary configurations con- 
sidered here, this factor is not the same across different spin 
configurations. However, because the £ — to = 2 represents 
by and large the most important contribution to the radiation 
emitted in the two channels and because the GW-emission 
from binaries with spins aligned/antialigned with the orbital 
angular momentum has been computed in a number of related 
works |f39l |40l |4Tl [42l [431 |43l |46l [471 [481 1491 [501 |5T1 [32l [53l 
|54| . the results found here allow us to simply extend all of the 
phenomenology reported so far for the £ — 2, m = 2 GW- 
emission from the above cited works also to the EM channel. 

Unfortunately, the tight correlation found in the amplitude 
evolution of the lowest-order mode disappears for higher- 
order modes. This is reported in the bottom panels of Fig. [8] 
which indicate that while the phase amplitude remains the 
same (cf. the bottom right panel), the evolution of the am- 
plitude in the two channels does not differ only by a simple 
constant scaling factor (cf. the bottom left panel). A sim- 
ilar behaviour is found for lower-order modes such as the 
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TABLE II: Relative emitted energies in EM waves and GWs 
(E^ /M, /M, respectively), and emitted angular momentum 
in GWs (J™J /M 2 ), for the magnetic field B M = 10 -4 . 



model 


EM / 




Jot/M 2 


S-6 


1.562E - 7 


0.0243 


0.216 


ro 


2.040E - 7 


0.0357 


0.213 


so 


2.055E - 7 


0.0354 


0.243 


S6 


3.412E - 7 


0.0590 


0.380 



I — 4, to — 2 one but is not reported here for compact- 
ness. Overall, these results suggest that although the main 
(and lowest-order) contribution to the EM emission does in- 
deed come as a result of the dragging of the EM fields by 
the orbital motion of the binary, additional contributions arise 
from higher-order modes which are not directly related to the 
orbital motions. These are likely to be the result of the com- 
plex interactions among the EM fields, discussed in Fig.[3]and 
whose investigation, although of great interest, goes beyond 
the scope of this paper. 

Another interesting quantity to consider in our analysis is 
the energy carried away from the systems in the two emis- 
sions, which can be computed by using eqs. (f27])-(|2"8"]>, where 
we have taken into account modes up to I = 4. Despite 
the differences between the EM and gravitational waves dis- 
cussed already, when looking at the emission in the lowest- 
order modes that can be associated to the different multipolar 
decomposition of the two emissions (cf. Fig. |7j, we find that 
the overall energy fluxes are extremely similar and differ es- 
sentially only by a constant (but large) factor. This is shown in 
Fig.[9]which reports both the GW (left panel) and the EM radi- 
ated energy fluxes (right panel) when integrated over a sphere 
located at r — 100 M for all the binary sequences considered 
here. Once again, the fact that F EM basically mimics F GW , 
underlines that the emission in the EM channel is intimately 
tied to the emission in GW, so that the observation of one of 
the two would lead to interesting information also about the 
other one. As a final comment it is worth noting that although 
the energy fluxes from the binaries Sq and r show a different 
evolution, the total emitted energy, namely the area under the 
curves, is extremely similar and is reported in Table [II] This 
provides yet an additional confirmation of the results already 
presented in refs 11381 |4T1 l42l l43l [54 1 for binaries with aligned 
spins and yields further support to the conjecture that when 
the initial spin vectors are equal and opposite and the masses 
are equal, the overall dynamics of the binary is the same that 
of the corresponding nonspinning binary. 



VI. ASTROPHYSICAL DETECTABILITY 

As discussed in the previous sections, the EM and GW ra- 
diation are tightly coupled and evolve on exactly the same 
timescales and with the same spectral distribution in fre- 
quency. The rates of loss of energy and angular momentum, 
however, are very different. This is summarized in Table [II] 
which reports the total energy radiated during the inspiral and 



merger in either EM waves or GWs (i.e., £?ad /m, E 1 ^ /M) 
and the angular momenta radiated in GWs (i.e., J^/M 2 ). 
From the values obtained, two interesting observations can be 
made. The first one is that the radiated EM energy is higher 
for binaries which lead to a more highly spinning final black 
hole. This is a consequence of these binaries merging with 
increasingly tighter orbits and at higher frequencies, which 
leads to stronger EM and GW fluxes. The second one has al- 
ready been mentioned in the previous Section and reflects the 
fact that the binaries ro and sq lead to the same energy emis- 
sion (and to the same final black-hole spin ATI l42l |43l [54)) 
despite the sq binary has black holes with non-zero individual 
spins. 

Note also that, in contrast with the losses in the GW emis- 
sion, those in the EM one do not depend just on the masses and 
initial spins of the black holes but also on the strength of the 
initial magnetic field. This dependence must naturally scale 
quadratically with the magnetic field, so that we can write 
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(29) 



= 1.43* 10-*, ,30, 
where we have used the following relation 



B [G] = 8.36 x 10 



19 



Me 
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B [geom. units] . (31) 



to convert a magnetic field in geometric units 
(B [geom. units]) into a magnetic field expressed in 
Gauss (B [G]). 

As discussed before, the EM emission is closely related via 
simple scaling factors to the GW one and whose efficiency has 
been discussed in detail in Sect. VB of ref. [38]. In particu- 
lar, it was shown there that the radiated GW energy depends 
quadratically on the total dimensionless spin (see eq. (24) 
in 11381 ) and the corresponding coefficients pi were presented 
in eq. (25) in the same reference. Hence, at least in the case of 
equal-mass binaries, it is trivial to express k\{a\, a 2 , Mi, M 2 ) 
in terms of the suitably rescaled coefficients pi in [38]. Here, 
however, because we are interested in much simpler order-of- 
magnitude estimates, we will neglect the dependence of ki on 
the spins and simply assume that fc x ~ 10~ 7 , so that 
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(32) 



where we have considered a total black hole mass of 10 s M Q 
and a magnetic field of 10 4 G as representative of the one pos- 
sibly produced at the inner edges of the circumbinary disc ll55ll 
(see ll56l for a recent discussion on the strength of magnetic 
fields in active galactic nuclei (AGN)). 

It should be noted that only when an extremely strong mag- 
netic field of ~ 10 11 G is considered, does the EM efficiency 
become as large as E*™r/M ~ 10 _1 and thus comparable 
with the GW one. For more realistic magnetic fields, how- 
ever, and assuming for simplicity that E r J^/M ~ 10~ 2 for 
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all possible spins, the ratio of the two losses is 
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That is, for a realist value of the initial magnetic field, the GW 
emission is 13 orders of magnitude more efficient than the EM 
one. More importantly, however, the frequency of variation of 
the EM fields is of the order 
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(34) 



and therefore much lower than what is accessible via astro- 
nomical radio observations, which are lower-banded to fre- 
quencies of the order of ~ 30 MHz. As a result, it is very 
unlikely that a direct observation of the induced EM emission 
would be possible even from this simplified scenario. 

Nevertheless, in the spirit of assessing whether this large 
release of EM radiation can lead to indirect observations of an 
EM counterpart, it is useful to compare with the typical 
luminosity of an AGN. To fix the ideas let us consider again 
a black hole of mass M = 10 s M Q ~ 10 41 g ~ 10 61 erg, so 
that the luminosity in EM waves for Bq = 10 4 G will be 



rad 



T — ^EM 



IQ^(JL; \ ergs- 1 



Vl0 4 G 

R f B 
10 8 ' 



Vl0 4 G 



L 



i 



10 \Wg 1 w 



(35) 



where we have assumed a timescale r ~ 10 M ~ 10 s ~ 
1 d and where L & , l^Edd are the total luminosity of the Sun 
and the Eddington luminosity LEdd = 3.3 x 10 4 (M/M & )Lq, 
respectively. While this is a rather small luminosity (distant 
quasars are visible with much larger luminosities of the order 
10 47 erg s _1 ), it is comparable with the luminosity of nearby 
AGNs and that is of the order of 10 41 erg s _1 . More impor- 
tant, however, is the comparison between the EM emitted by 
the merging binary and the one coming from the accretion 
disc. Using ( |3~5| ) it is straightforward to deduce that the binary 
EM luminosity is comparable with that of an AGN accreting 
at 10~ 4 the Eddington rate. Hence, unless the accretion rate is 
rather small (namely, much smaller than 10~ 4 the Eddington 
rate with the extreme case being the non-accreting scenario) 
the EM emission from the binary would be not only restricted 
to very low-frequencies but also just a small fraction of the to- 
tal luminosity. Under these conditions it is unlikely that such 
emission could have an observable impact on the overall lu- 
minosity of the accreting system. 

As a final consideration it is useful to estimate whether the 
inspiralling binary could nevertheless imprint a detectable ef- 
fect on the disc via the perturbations in the magnetic field it 
can produce. To assess whether this is the case we first com- 
pare the frequency / B with the typical plasma frequency 
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where n e is the electron number density, or with the electron 
cyclotron f c frequency 



fc 



lu c eB 
2tt 2nm P c 



10 



10 



( B 



Vl0 4 G 



Hz. 



(37) 



Clearly, the magnetic field varies with a frequency / B that is is 
between 14 and 18 orders of magnitude smaller and hence that 
the electrons and protons in the disc are always able to "ad- 
just" themselves to the changes in the magnetic fields, which 
are extremely slow when compared with the typical timescales 
in the plasma. Stated differently, the EM radiation produced 
by the inspiral cannot penetrate the disc and will be effectively 
reflected over a skin depth of A = c/uo e ~8x 10~ 6 cm. 

Finally, we consider whether the perturbed magnetic mag- 
netic field can have impact on the transport of angular mo- 
mentum in the disc and hence modify its accretion rate in a 
detectable way. It is worth remarking, in fact, that there is con- 
siderably large EM energy flux reaching the accretion disc and 
that is 



10 n (B/10 4 G) 2 erg s" 1 cnr 2 , 
where r; n ~ 10 z r g is the inner radius of the disc and 
r g ~ 10 15 cm is the gravitational radius for a black hole of 
10 8 M®. A crude way to estimate the perturbation on the 
disc is by considering the ratio between the viscous transport 
timescale Ty and the magnetic transport timescale induced by 
the oscillating magnetic field, tb- Should this ratio be of the 
order of unity (or larger), then the magnetic-field perturbation 
may be transmitted to the disc in the form of Alfven waves. 
In practice we estimate this by considering the (inverse) ra- 
tio between the viscous and magnetic torques, with the first 
one being expressed in terms of the average pressure p and 
sound speed c s as f^y ~ ap ~ oipc 2 s and the second one 
as f^ B ~ r(5B*B z a/(87r) ~ r/3S 2 /(87r); here a is the 
standard alpha-disc viscosity parameter and is a measure 
of the perturbation induced in the background magnetic field 
(i.e., SB^ ~ (3B , B z ~ B ). We therefore obtain 
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where 10 



is the typical length scale over which mag- 



netic torques could operate. Assuming now a ~ 0.1 — 0.01, 
(3 ~ 10~ 2 and c s 0.1 — 0.01c as reference numbers, the 
rough estimate ( [38) suggests that it is indeed possible that 
tv > tb and hence that the perturbations in the magnetic 
field, albeit small and rather slow, can induce a change in the 
viscous torque and hence induce a change in the accretion rate 
if the latter is sufficiently stable. Determining more precisely 
whether this modulation in the magnetic field can effectively 
leave an imprint on the accretion flow would require a more 
accurate modeling of the accretion disc and is clearly beyond 
the scope of this simple estimate. It is however interesting that 
this possibility is not obviously excluded. 

In summary, the analysis carried out in this Section shows 
that it is highly unlikely that the EM emission associated with 
the scenario considered in this paper can be detected directly 
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and simultaneously with the GW one. This is essentially be- 
cause the EM is too inefficient for realistic values of the mag- 
netic fields and because it operates at frequencies which are 
well outside the ones accessible to astronomical radio obser- 
vations. However, if the accretion rate of the circumbinary 
disc is sufficiently stable over the timescale of the final inspiral 
and merger of the black-hole binary, then it may be possible 
that the EM emission will be observable indirectly as it will 
alter the accretion rate through the magnetic torques exerted 
by the distorted magnetic field lines. A firmer conclusion of 
whether this can actually happen in practice will inevitably 
have to rely on a more realistic description of the accretion 
process. 

As a final comment we stress that our analysis and discus- 
sions have not included the role of gas or plasmas around the 
black hole(s) nor have we considered resistive scenarios. Both 
of these ingredients, when coupled to the EM fields behavior 
described here, could induce powerful emissions by acceler- 
ating charged particles via the strong fields produced (e.g., in 
a manner similar to the Blandford-Znajek mechanism |[T6l ) or 
by affecting the gas/plasmas dynamics or via the reconnec- 
tion of the complex EM fields produced during the inspiral 
and merger. Future work in these directions is needed in or- 
der to shed light on these possibilities and assess their realistic 
impact as EM counterparts to the GW emission. 

VII. CONCLUDING REMARKS 

We have analyzed the phenomenology that accompanies 
the inspiral and merger of black-hole binaries in a uniform 
magnetic field which is assumed to be anchored to a distant 
circumbinary disc. Our attention has been concentrated on 
binaries with equal masses and equal spins which are either 
aligned or antialigned with the orbital angular momentum; in 
the case of supermassive black holes, these configurations are 
indeed expected to be the most common ones [19 20 1. Fur- 
thermore, this choice allows us to disentangle possible pre- 
cession effects and concentrate on the EM fields dynamics as 
affected by the orbital motion of the binary. Overall, the sim- 
ulations reveal several interesting aspects in the problem: 

• The orbital motion of the black holes distorts the essen- 
tially uniform magnetic fields around the black holes 
and induces a quadrupolar electric field analogous to 
the one produced by the Hall effect for two conduc- 
tors rotating in a uniform magnetic field. In addition, 
both electric and magnetic fields lines are dragged by 
the orbital dynamics of the binary. As a result, a time 
variability is induced in the EM fields, which is clearly 
correlated with the orbital behavior and ultimately with 
the GW-emission. The EM fields become, therefore, 
faithful tracers of the spacetime evolution. 

• As a result of the binary inspiral and merger, a net 
flux of electromagnetic energy is induced which, for 
the I — 2, to ± 2 modes is intimately tied, via a con- 
stant scaling factor in amplitude, to the gravitational en- 
ergy released in GWs. This specular behaviour in the 



amplitude evolution disappears for higher-order modes, 
even though the phase evolution remains the same for 
all modes. 

• Because the tight correlation between the EM and the 
GW-emission has been found for all of the cases consid- 
ered here, we expect it to extend to all possible binary 
configurations as long as the EM fields are playing the 
role of "test fields". Hence, the modelling of the GW 
emission does in practice provide information also on 
the EM one within the scenario considered here. 

• Although the global large-scale structure of the EM 
fields is dictated by the orbital motion, the individual 
spins of the black holes further distort the EM field lines 
in their vicinities. These small-scale fields may lead to 
interesting dynamics and to the extraction of energy via 
acceleration of particles along open magnetic field lines 
or via magnetic reconnection. 

• The energy emission in EM waves scales quadrati- 
cally with the total spin and is given by E T ^~/M ~ 

1CT 15 (Af/10 8 A/©) 2 (B/10 4 G) 2 , thus being 13 or- 
ders of magnitude smaller than the gravitational energy 
for realistic magnetic fields. This EM emission is at fre- 
quencies of ~ 1CT 4 (1O 8 M /Af) Hz, which are well 
outside those accessible to astronomical radio observa- 
tions. As a result, it is unlikely that the EM emission 
discussed here can be detected directly and simultane- 
ously with the GW one. 

• Processes driven by the changes in the EM fields could 
however yield observable events. In particular we ar- 
gue that if the accretion rate of the circumbinary disc 
is small and sufficiently stable over the timescale of the 
final inspiral, then the EM emission may be observable 
indirectly as it will alter the accretion rate through the 
magnetic torques exerted by the distorted magnetic field 
lines. 

All of these results indicate that the interplay of strong grav- 
itational and EM fields represents a fertile ground for the de- 
velopment of interesting phenomena. Although our analysis 
is incomplete as the effects on plasmas are not taken into ac- 
count, we believe that the main properties of the EM dynamics 
described above should hold as long as the energy in the black 
holes dominates the energy budget. A more precise estimate 
of the possible emissions and of the observational signatures 
calls for further studies which would necessarily have to in- 
clude additional physics. This work, however, together with 
those in refs. El [EJ [13] [T4J [151 [121 [lUl , constitute interesting 
first steps in this direction. 



Acknowledgments 

It is a pleasure to thank Kyriaki Dionysopoulou, Olindo Zan- 
otti, Chris Thompson, Avery Broderick, Steve Liebling and 
David Neilsen for useful discussions and comments. DP has 



15 



been supported as VESF fellows of the European Gravita- 
tional Observatory (EGO). Additional support comes from 
the DFG grant SFB/Transregio 7 "Gravitational Wave As- 
tronomy", NSF grant PHY-0803629 and NSERC through a 
Discovery Grant. Research at Perimeter Institute is sup- 



ported through Industry Canada and by the Province of On- 
tario through the Ministry of Research & Innovation. The 
computations were performed at the AEI, the LONI network 
(www . loni . org), LRZ Munich, and Teragrid. 



[1] B. S. Sathyaprakash and B. F. Schutz, Living Rev. Rel. 12, 2 

(2009), 0903.0338. 
[2] J. S. Bloom et al. (2009), 0902.1527. 

[3] P. J. Armitage and P. Natarajan, Ap. J. Letters 567, L9 (2002). 
[4] F. K. Liu, X.-B. Wu, and S. L. Cao, MNRAS 340, 411 (2003). 
[5] M. Milosavljevic and E. S. Phinney, Astrophys. J. 622, L93 

(2005), astro-ph/04 10343. 
[6] S. M. O'Neill, M. C. Miller, T. Bogdanovic, C. S. Reynolds, 

and J. Schnittman, ArXiv e-prints (2008), arXiv:08 12.4874. 
[7] L. R. Corrales, Z. Haiman, and A. MacFadyen (2009), 

0910.0014. 

[8] J. D. Schnittman and J. H. Krolik, Astrophys. J. 684, 835 

(2008), 0802.3556. 
[9] Z. Lippai, Z. Frei, and Z. Haiman, Astrophys. J. Lett. 676, L5 

(2008), 0801.0739. 
[10] G. A. Shields and E. W. Bonning, Astrophys. J. 682, 758 

(2008), 0802.3873. 
[11] M. Megevand, M. Anderson, J. Frank, E. W. Hirschmann, 

L. Lehner, S. L. Liebling, P. M. Motl, and D. Neilsen, Phys. 

Rev. D. 80, 024012 (2009), 0905.3390. 
[12] A. I. MacFadyen and M. Milosavljevic, Astrophys. J. 672, 83 

(2008) . 

[13] P. Chang, L. E. Strubbe, K. Menou, and E. Quataert (2009), 
0906.0825. 

[14] J. R. van Meter et al. (2009), 0908.0023. 

[15] T. Bode, R. Haas, T. Bogdanovic, P. Laguna, and D. Shoemaker 

(2009) , 0912.0087. 

[16] R. D. Blandford and R. L. Znajek, Mon. Not. Roy. Astron. Soc. 
179, 433 (1977). 

[17] C. Palenzuela, M. Anderson, L. Lehner, S. L. Liebling, and 

D. Neilsen, Phys. Rev. Lett. 103, 081101 (2009), 0905.1121. 
[18] C. Palenzuela, L. Lehner, and S. Yoshida (2009), 091 1.3889. 
[19] T. Bogdanovic, C. S. Reynolds, and M. C. Miller, ApJ 661, 

L147 (2007), arXiv:astro-ph/0703054. 
[20] M. Dotti, M. Volonteri, A. Perego, M. Colpi, M. Ruszkowski, 

and F. Haardt, MNRAS pp. 1795-+ (2009), 0910.5729. 
[21] C. Palenzuela, L. Lehner, O. Reula, and L. Rezzolla, 

arXiv:0810.1838 (2008), 0810.1838. 
[22] S. S. Komissarov, MNRAS 382, 995 (2007). 
[23] D. Pollney, C. Reisswig, E. Schnetter, N. Dorband, and P. Di- 

ener (2009), 0910.3803. 
[24] L. Baiotti, B. Giacomazzo, and L. Rezzolla, Phys. Rev. D 78, 

084033 (2008), ArXiv e-prints 0804.0594. 
[25] http://www.cactuscode.org. 

[26] E. Schnetter, S. H. Hawley, and I. Hawke, Class. Quantum Grav. 

21, 1465 (2004), gr-qc/03 10042. 
[27] D. Pollney, C. Reisswig, N. Dorband, E. Schnetter, and P. Di- 

ener (2009), 0910.3656. 
[28] L. Gunnarsen, H. Shinkai, and K. Maeda, Class. Quantum Grav. 

12, 133 (1995), gr-qc/9406003. 
[29] A. M. Abrahams, L. Rezzolla, M. E. Rupright, and et al., Phys. 

Rev. Lett. 80, 1812 (1998), gr-qc/9709082. 
[30] M. E. Rupright, A. M. Abrahams, and L. Rezzolla, Phys. Rev. 

D 58, 044005 (1998). 



[31] L. Rezzolla, A. M. Abrahams, R. A. Matzner, M. E. Rupright, 
and S. L. Shapiro, Phys. Rev. D 59, 064001 (1999), gr- 
qc/9807047. 

[32] A. Nagar and L. Rezzolla, Class. Quantum Grav. 22, R167 

(2005), erratum-ibid. 23, 4297, (2006). 
[33] L. Lehner and O. M. Moreschi, Phys. Rev. D76, 124040 (2007), 

0706.1319. 

[34] R. M. Wald, Physical Review D 10, 1680 (1974). 

[35] R. Crowley, D. Macdonald, R. Price, I. Redmount, Suen, 

K. Thorne, and X.-H. Zhang, Black Holes: The Membrane 

Paradigm (Yale University Press, 1986). 
[36] M. Ansorg, B. Briigmann, and W. Tichy, Phys. Rev. D 70, 

06401 1 (2004), gr-qc/0404056. 
[37] L. Rezzolla et al, Astrophys. J 679, 1422 (2008), 0708.3999. 
[38] C. Reisswig et al. (2009), 0907.0462. 

[39] M. Koppitz, D. Pollney, C. Reisswig, L. Rezzolla, J. Thornburg, 
P. Diener, and E. Schnetter, Phys. Rev. Lett. 99, 041 102 (2007), 
gr-qc/0701163. 

[40] D. Pollney, C. Reisswig, L. Rezzolla, B. Szilagyi, M. Ansorg, 
B. Deris, P. Diener, E. N. Dorband, M. Koppitz, A. Nagar, et al., 
Phys. Rev. D76, 124002 (2007), arXiv:0707.2559. 

[41] L. Rezzolla, E. N. Dorband, C. Reisswig, P. Diener, D. Poll- 
ney, E. Schnetter, and B. Szilagyi, ApJ 679, 1422 (2008), 
0708.3999. 

[42] L. Rezzolla, P. Diener, E. N. Dorband, D. Pollney, C. Reis- 
swig, E. Schnetter, and J. Seiler, ApJ 674, L29 (2008), 
arXiv:0710.3345. 

[43] L. Rezzolla, E. Barausse, E. N. Dorband, D. Pollney, C. Reis- 
swig, J. Seiler, and S. Husa, Physical Review D 78, 044002 

(2008) , 0712.3541. 

[44] S. Husa, M. Hannam, J. A. Gonzalez, U. Sperhake, and 

B. Bruegmann, Phys. Rev. D77, 044037 (2008), 0706.0904. 
[45] M. Campanelli, C. O. Lousto, and Y. Zlochower, Phys. Rev. D 

74, 041501 (2006), gr-qc/0604012. 
[46] M. Campanelli, C. O. Lousto, and Y. Zlochower, Phys. Rev. D 

74, 084023 (2006). 
[47] M. Campanelli, C. O. Lousto, and Y. Zlochower, Phys. Rev. D 

73, 061501 (2006), gr-qc/0601091. 
[48] A. Buonanno, L. E. Kidder, and L. Lehner, Physical Review D 

77, 026004 (2008), 0709.3839. 
[49] F. Herrmann, I. Hinder, D. M. Shoemaker, P. Laguna, and R. A. 

Matzner, Physical Review D 76, 084032 (2007), 0706.2541. 
[50] F. Herrmann, I. Hinder, D. Shoemaker, P. Laguna, and R. A. 

Matzner, Astrophys. J. 661, 430 (2007), arXiv:gr-qc/0701 143. 
[51] L. Boyle, M. Kesden, and S. Nissanke, Physical Review Letters 

100, 151 101 (2008), arXiv:0709.0299. 
[52] L. Boyle and M. Kesden, Phys. Rev. D78, 024017 (2008), 

0712.2819. 

[53] P. Marronetti, W. Tichy, B. Bruegmann, J. Gonzalez, and 

U. Sperhake, Phys. Rev. D77, 064010 (2008), 0709.2160. 
[54] E. Barausse and L. Rezzolla, Astrophys. J. Lett. 704, L40 

(2009) , 0904.2577. 

[55] M. Milosavljec and E. S. Phinney, Astrophys. J. 622, L93 
(2005). 



16 



[56] N. A. Silant'ev, M. Y. Piotrovich, Y. N. Gnedin, and T. M. 
Natsvlishvili, ArXiv e-prints (2009), 0908.2725. 



